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A s=^ antiferromagnetic spin chain is equivalent to the two- flavor massless Schwinger model 
in an uniform background charge density in the strong coupling. The gapless mode of the 
spin chain is represented by a massless boson of the Schwinger model. In a two-leg spin 
ladder system the massless boson aquires a finite mass due to inter-chain interactions. The 
gap energy is found to be about .25fc \J'\ when the inter-chain Heisenberg coupling J' is small 
compared with the intra-chain Heisenberg coupling, fc is a constant of 0(1). It is also shown 
that a cyclically symmetric Af^-leg ladder system is gapless or gapful for an odd or even A'^^, 
respectively. 



A s = i spin chain with antiferromagnetic nearest-neighbour Heisenberg couphngs is 
exactly solved by the Bethe ansatz.[|I| It has a gapless excitation. A two-leg spin ladder 
consists of two spin chains coupled each other. Experimentally a two-leg spin ladder system 
has no gapless excitation. |Hll The gapless mode of spin chains becomes gapful. In this 
paper we give, without resorting to numerical evaluation, a deductive microscopic argument 
which shows why and how it happens. 

Spin ladder systems are not exactly solvable. Various approximation methods have been 
employed in the literature. 0, ||, ||, ^ |ll| We first show that a s = ^ spin chain is 
equivalent to the two-fiavor massless Schwinger model in the strong coupling in an uniform 
background charge density. The two-flavor Schwinger model has a massless boson excitation, 
which corresponds to the gapless excitation in the Bethe ansatz. A spin ladder system is 
described as two sets of two-flavor Schwinger models which interact with each other by 
four-fermi interactions. 

An antiferromagnetic spin chain is described by 

HchMS) =jJ2Sn- Sn+1 (J > 0) (1) 

wheras a two-leg spin ladder is described by 

^^laddor(S', T) = -ff chain (S') + i?chain(T) + Hrung{S , T) 
Hi-ung{S,T) = j'^SVi'^n- (2) 

Consider first a, s — anti- ferromagnetic spin chain -ffchain(»5). We express the spin operator 
in terms of electron operators by Sn — cl^^(jc„. With the aid of the Fierz transformation 
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we have 45„ • Sn+i = -{4iC„+i, cl+iC„} + 1 - (4,c„ - l)(4j+iC„+i - 1). The last term 
drops as the half-filhng condition cJjC„ = 1 is satisfied for a spin chain. The first term is 
linearized by introducing an auxiriary field, or by the Hubbard-Stratonovich transformation. 
The Hamiltonian iJchain is equivalent to the Lagrangian 

-^chain = ^ I — A„(cJjC„ — 1) — -^{UnUn — C/„4iC„^_i — C^,tcli+lC„) | . (3) 

U„ is a link variable, defined on the link connecting sites n and n + 1. A„ is a Lagrange 
multiplier enforcing the half-filling condition at each site. The transformation is valid for 
J > 0. The Lagrangian i^pj^ has local U{1) gauge invariance. 

We consider a periodic chain of N sites: Sjy+i = Si. The mean field energy is evaluated, 
supposing \Un\ = U, to be -Bmcan = J{^NU^ - 2Ucot{-K/N) + \N}. For large N it has a 
sharp minimum a,i U ~ 2/-K. Radial fiuctuations of t/„'s are suppressed, though quantum 
fiuctuations of the phase of C/„'s cannot be neglected. We write 

f/„ = -e^^^" (4) 

TT 

where I is the lattice spacing. We need to incorporate quantum fiuctuations of A„ and An to 
all orders. With (^) substituted the Lagrangian (||) becomes that of lattice electrodynamics. 

To make this point clearer, we take the continuum limit. For an anti-ferromagnetic spin 
chain, two sites form one block. The even-odd site index becomes an internal (spin) degree 
of the Dirac field in the continuum limit. The correspondence is given by 

'\l)'f'\x) = -1= — C2s~iM at odd site 



{x) — — 1=^ C2s,a at even site 



2£ 
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where x corresponds to (2s— 1)1' and 2sl. With the given normalization {V'j"'' '^'k^ (j/)^} = 
5°''^5jkdL{x — y) in the continuum limit, where 5]^{x) is the periodic delta function with the 
period L = N£. The phase factors in (||) refiect the fermi momentum kp = i^vr at the half 
filling. 

The term J2n '4i^n+i + ^-C- becomes 2i£J2a I dx(->p[''^^dxip2^^ + ip2^^'^ d^'tpi'^) ■ Hence in 
the continuum limit the original spin Hamiltonian (|l|) is transformed to a system with the 
Lagrangian density 

2 

>Ce\ai„[^M:^] = -^^M. +EcV^^"^7^(*?i5^ - lA,y^^ + . (6) 

a—1 

Here the Dirac matrices are 7" = (T3, 7^ = i(T2- The "light" velocity c is given by c = 2£J/ttTi. 
Xq = ct and (Aq, Ai) = {\,cA). Although the Maxwell term is absent in the £ ^ limit, 
it is generated at finite I. The coupling constant e must be expressed in terms of J and I. 
From the dimensional analysis = k'^J/t where fc is a constant of 0(1). Note that in the 
t Q limit with c kept fixed, diverges as 
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This is nothing but the two-flavor massless Schwinger model in the strong couphng in 
an uniform background charge density. The term representing the background charge 
arises from the half-fiUing condition. The system is neutral as a whole. 

Notice that the spin index a of original electrons becomes a flavor index in (|^), while 
the even-odd index j becomes a spin index of the Dirac fleld ij}^°''^{x). The two-flavor nature 
reflects the electron spin i. 

The correspondence of the spin chain model to QED2 has been noted in the literature, 
but the rigorous derivation has not been given before, fl^ In particular, the importance of 
the two-flavor nature has not been recognized. Mapping to SU{2) gauge theory also has 
been suggested. 

The two-flavor massless Schwinger model is exactly solvable, [m] Quantum fluctuations of 
aU flelds, ip'^'''^ and A,,, can be completely taken into account. With the periodic boundary 
condition, the model is two-flavor QED2 deflncd on a circle, which has been analysed in 
detail by the bosonization method. ||l5|, p^ , 

The bosonization formula for the left and right moving components of the Dirac flelds is 

V-il^,.!) = ^CJe±*{9±+2''P±(*±^)/-^>7Vo[e±*^"^±(*'^)] (a = 1,2) (7) 

where = e"SLi(P++P-) and CI = e'''Sb=i(P+"P-\ (0^) represents left (right) 
moving modes. Nq \ ] denotes the normal ordering in a basis of massless fields. The Hamil- 
tonian becomes jlTl 

a=l 

+ f'dx^{^^' + ^'' + ^'i>' + ^x' + x") ■ (8) 
Jq 2 ^c/ nhc 

The neutrality condition reads Qi + Q2 = L/l — N . Qw and Pw are the Wilson line 
phase e'^"' = exp \(i/hc) Jq dx Aij and its conjugate momentum. Qa = —pt + Pa and 
Qsa = Pa + Pa charge and axial charge of the a-th flavor, respectively, both of which 
take integer eigenvalues and commute with the Hamiltonian. $ = (0i + (f)2)/V2 and x — 
(01 — 4>2)/ where (f)a — + (j)"_ and / dxcf)a — 0. 

The $ field has a Schwinger mass ji where /i^ — 2e'^'h/Tic'. The excitation energy is 
li<? = 2kJ/'K ^ .637k J. The x field is massless, which corresponds to the gapless excitation 
in the spin chain and controls the behavior of correlation functions at large distances. The 
wave function for the zero mode part is written as 

I*) = E / '^P^ \pw,n,r)e-'^'^+^^"'P-'fipw,V + 7Tpw) 

n,r 

Pw\pw,n,r) = pw\pw,n,r) 

p±\pw,n,r) = {n + rda,i T \N)\pw,n,r) (9) 
where f{pw,v) uiust solve the Schrodinger equation 

K{pw,f)f{pw,(p) ^<^f{vw,^) 
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For the ground state f{pw, f) = const • g—^t^'^i^Pw/^^f^, 

In the Schwinger model there is a parameter characterizing states. The wave func- 
tion (^) corresponds to = 0. The 9 vacuum originates from the invariance under large 
gauge transformatins and the chiral anomaly in the continuum theory, jl^ In the lattice spin 
systems the lowest energy state with 6* = is expected to be singled out. 

Employing the bosonization formula, the critical exponent of the spin-spin correlation 
function ( S{2n)S{0) ) ^ n^'^ (n ^ 1, n <C N) is found to be 77 = 1, which agrees with the 
result from the Bethe ansatz. 

Now we consider a spin ladder system (|^). In the absence of the inter-chain rung in- 
teraction (J'=0) the system is equivalent to the two sets of two-flavor massless Schwinger 
models described by Cli^^^^[Afj,,'ip] -I- /Z^j^j^j^JA^, , V'] . With the aid of the correspondence (||), 
the inter-chain interaction i?rung in the continuum limit is written as 

J' N tJ' f / ~~ ~~ _~— _~~\ 

^™ng = — + — + + + (11) 

where every quantity in the expression is flavor singlet; ip^ip ~ 5^L=i ''/''■'''' etc. Notice 
that both charge density and scalar density operators appear in (O). The chiral symmetry 
is broken, which leads to mass generatioii. 

When expressed in terms of Tp± and ip±j ^rang contains many terms. The Hamiltonian 
is simplified in the large volume limit L ~ N£ — > 00. Define pa — ip^""^^ Ma = V'+''^V'-"\ 
and corresponding pa and Ma- Relevant terms in H^^^^ are 

-^rung ~ + H^b 



J'£ f 

H^a = — j dx {pi - P2)[pi - P2) 



J't 

Hsb = — I dx 



■ [{Ah - M2){mI - Ml) + (h.c.)} . (12) 



Terms of the form MaMf, are suppressed as fluctuations in Qa are small compared with the 
average N/2. 

Boson fields associated with ip and ip are denoted by (<&, x) and ($, x), respectively. We 
introduce a new orthonormal basis: $± = (<I> ± $)/V2 and x± — ix ^ x)/^- The first 
term in (|l^ ) is 



(Qi - Q2)(Qi - O2) + / dx—{{d.,x+f - id.X-f} ■ (13) 



It changes the propagation velocities of x± fields. 
It follows from (^ that 

M^M^ = e^^'^*('5"^'5'')^/-^e~*'^'"~''''^ — A^o[e^*^^^'^°~'^''^] • (14) 
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Note that No[el^^] = B{mcL/h)>^^ Z'^'' Nm[e"^^] where the reference mass in the normal 
ordering N[ ] is shifted from to m. B{0)=1 and B{z) ^ e^zj^-n for z > That 
is, if all fields become massive, ( [T4| ) is nonvanishing in the L oo limit. Otherwise 
vanishes. In passing, terms not included in are suppressed exponentially in the L —i- oo 
limit when x± fields aquire masses. 

There are fiuctuations in Qa- Write Qi^2—-^N zLQ and Qi^2—^N±Q. Important terms 
in f dxMaMl result when Q = ±Q. Since \Q\, \Q\ <^ N, we have in the large volume limit 



4 KAttH/ 

-/x<i>_/Xx+{e-^(«i-''2)^[g-«v^(*-.+x+)] ^_g-'(92-gi)^[g-«v^(<E._-x+)]| 



(15) 



Here we have defined $± = ($ ± <I)/V2 and x± = (x ± x)/V2. iV[e"*^(*-+x-)] denotes 
that the <&_ and X- fields are normal-ordered with respect to their masses and , 
respectively. 

has two major effects. It gives an additional potential in the zero mode sector: 

x{mx- K'^^i-'^i) +e-*(«^-«^)] ~ ^i^^[e"'^'^'-^'^ + e-'^'^^--'^'^] + (h.c.)J . (16) 
Secondly it gives additional masses to and x±- For small | J'| ^ J 



2 



I? - „ /p±' 91-91 \ 

/^L = £S^*-^x.(e^^^---^) (17) 

Here we have made use of ( e±'(?i-«i) > = ( e^'^'^-gs) ) and ( e±'(?i-«=^) ) = (e±'fe-«i)), 
which reflects the up-down symmetry of the original spin system and is justified shortly. 

The wave function of the ladder system is specified with f{pw, 'P'lPw, 0) as in (||). The 
rung interaction (|l^) gives an additional potential in the ip representation, e**^ and e*'^ give 
rise to q^'p-^^pw ^nd q-^v-^^pw ^ respectively. / satisfies 



K-ung = cos(v3 - (^) - A*x+ cos(v3 + <^)| C0S7r(pvF -Pw) • (18) 
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For large L the potential term dominates in Eq. ([l8|). The ground state wave function has 
a sharp peak at the minimum of the potential. For J' > (J' < 0), the minimum occurs at 
Pw = Pw = and ip = —(p = ±^7r {(f = (p = i^Tr) so that 

^g±«(9„-g„)^ ^ _^g±z(gi-q2)) ^ _(g±»fe-9i)) = for { ^' > (19) 

The masses are determined by (0) and ([l9|): 
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Vl - 2k2 '^^+ VI - 2k 

47r ?ic " 8 J "-^^^ J • ^ 
The expression is valid for small k. The excitation energy, a spin gap, is 

Aspin = Mx±c' ~ «;Aic2 = -— I J'l = 0.25 fc I J'l . (21) 

47r 

The ratio of Agpin to /zc^ is k. The gapless mode becomes gapful. The spin gap is determined 
by I J' I, generated irrespective of the sign of J'. The energy density is lowered: 

A£ = . (22) 

We have shown that the rung interaction breaks the chiral symmetry of spin chain 
systems, and generates a spin gap. 

In the literature the spin gap has been determined by various numerical methods for 
varying J'/J.Q In particular, Greven et al. obtained Agpin = AIJ' for small J' / J and .50 J' 
for J' — J, which is consistent with our prediction (pT|). 

It has been well known that spin chain systems are mapped to non-linear sigma models. ||l8|| 
Sierra has applied this mapping to iV^-leg ladder systems of spin S, and has shown that the 
spectrum is gapful or gapless for an integer or half-odd- integer SNi, respectively. The 
mapping to sigma models is valid for large SN^ 3> 1, while our method of mapping to the 
Schwinger model works for S = ^. 

The method of bosonization has been employed in the spin ladder problem. Schulz, in 
analysing a spin S chain, expressed 5 as a sum of 25* spin ^ vectors, thereby transforming 
the spin chain to a special kind of a spin i ladder system. With the aid of bosonization and 
renormalization group analysis he concluded that the spectrum is gapless for a half-odd- 
integer S.^ 

More recently a 2-leg s = 5 ladder system has been analysed by bosonization by Shelton 
et al. and by Kishine and Fukuyama. [pT| They have obtained a similar Hamiltonian to ours, 
but could not determined the gap. Our bosonization formula (|^) is a rigorous operator 
identity with no ambiguity in normalization, with which the Hamiltonian is transformed in 
the bosonized form. The correct treatment of the normal ordering is crucial in dealing with 
the mass (gap) generation. Not only the light modes (x±) but also the heavy modes (<&±) 
and zero modes {Q,qa) play an important role, which has been dismissed in ref. 11. 



6 



Our argument can be generalized to A^f-leg s — ^ ladder systems. Inter-chain interactions 
are given by i/iung = 12{ij) J'ij Sn^ where i and j are chain indices and (ij) labels 

rung pairs. J-^ = 2J for all (ij) in Schulz' model in ref. 10. 

Let us consider a cyclically symmetric antiferromagnetic ladder system in which non- 
vanishing Jlj's are J^i^i — J' > where J^^ = ^. Among boson fields $i's or 
Xi's, the singlet combination is denoted by $+ or x+- Other combinations of (f>'s or x's are 
degenerate. There are four masses to be determined: and /ix±- A'^i ^ M for small | J'|. 
The issue is whether or not all x fields become massive. The crucial part is the mass of x+- 

Repeating the above argument, one finds that the part of the rung potential K-ung in 
(|l8|), fj,^_ cos{ip — ip) — fi^^ cos{tp + (f), is replaced by 

i=l i=l 

where v^Aff+i = Vi- If Mx- ~ K-ung = and no correction arises to or /^^i- This 
solution has a higher energy density than the non-trivial solution so that ^ 0. From 
the symmetry Kung is minimized at cos((^i — ^i+i) = /- = Ij ■ ■ ■ : -^f). This implies that 
ipj = if + [j — l)ri and rj = 2p'K /N^ or 2p7r/(iVf — 2) where p is an integer. 

Suppose /ijj.^ ^ 0. Then cos((^i = /+ (i = 1, • ■ ■ , A'f). This leads to an additional 

condition that rj ^ n. All of these conditions are satisfied for an even Ni. The potential 
is minimized at f2p+i = i^Ti" and (p2p — Tj'"'- For an odd Ng the conditions cannot be 
satisfied. 

If ^-^^ = 0, 1] need not be tt. This gives a solution for an odd Nf. For an even iV^, 
this solution yields a higher energy density than the solution with ^ above. To 
summarize, the spectrum is gapless for an odd Ng, but is gapful for an even Ng. The 
interaction is frustrated in the rung direction for an odd Ng. The argument here is similar 
to Schulz' in ref. 10. 

In the experimental samples]^ J' ^ J so that k = 0(1). For instance, in SrCu203 (2-leg 
ladder), J ^ J' ^ 1300K and A^pin - 420K. The formula Q need to be improved by 
taking account of effects of nonlinear terms in ([l5|). Further it is observed that spin ladder 
systems with three legs are gapless. [The experimental sample is not cyclically symmetric: 
J'12 — J23 ^ J but = 0.] For this the large value of k is important, as our analysis 
indicates that a gap is generated so long as k is sufficiently small. It has been also reported 
that the spin gap is not affected by nonmagnetic impurities. We will come back to those 
points in separate publications. 
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